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ABSTRACT 


Various  facets  of  thermal  testing  of  space  vehicles  in  cold  vacuum 
chambers,  with  and  without  solar  simulation,  are  discussed.  Emphasis 
is  on  thermal  modeling  of  spacecraft.  Thermal  modeling  is  discussed 
for  the  case  of  simultaneous  preservation  of  temperature,  material  and 
surface  properties.  The  following  specific  items  are  treated:  (1)  a 
criterion  for  the  application  of  isothermal  modeling  is  developed  and 
illustrated  with  experimental  data,  (2)  rules  for  nonisothermal  modeling 
are  discussed  and  illustrated  with  experimental  results,  (3)  the  effect  of 
internal  convection  is  calculated  and  modeling  rules  are  developed, 

(4)  the  effect  of  internal  radiation  is  discussed  and  the  relative  importance 
illustrated  theoretically,  and  <5)  a  criterion  is  developed  for  the  degree  of 
intensity  nonuniformity  for  a  solar  simulator  in  terms  of  desired  tempera¬ 
ture  accuracy. 
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NOMENCLATURE 

Dimensionless  amplitude  of  solar  simulator  nonuniformity 
Surface  area,  m2 
Reference  area,  m2 

ceMe 

Ratio  of  heat  capacities,  .. 

csMa 

Specific  heat,  Joule/kg ~°K. 

Material  thickness,  m 
Nondimensional  temperature  error 
Radiation  function  in  Eq.  (19) 

Heat  transfer  coefficient,  watt/m2-°K 
Heat  transfer  coefficient,  watt/m2  -  °K 
Reference  heat  transfer  coefficient,  watt/m2  -  °K 
Nondimensional  heat  transfer  coefficient 
Scaling  ratio 

Reference  length,  linear  dimension  of  body,  m 
Half  wavelength  of  solar  simulator  nonuniformity,  m 
Mass,  kg 
Integer 

^  (Te  -  T0)s 
T„ 

te-t0 

Heat  flux,  watt/m2 

q 

Nondimensional  heat  flux,  — T  - 
Radius,  m 

Internal  heat  generation,  watt/m3 

Temperature,  °K 

Time,  sec 

Volume,  ml 

Reference  volume,  m3 
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Cartesian  coordinates,  m 
Absorptivity 
Nondimensional  area 
Nonuniformity  parameter,  Eq.  (21) 

Conduction  parameter,  Eq.  (19) 

Ratio  of  emissivities,  tj/fB 
Surface  emissivity 
Nondimensional  temperature 
Thermal  conductivity,  watt/m  -  :K 
Cos  i/ 

Polar  angle  for  sphere 

Nondimensional  coordinates  in  x-,  y-,  and  z-direction 
Density,  kg/m1 

Stefan-Boltzmann  constant,  watt/m1- (®K)4 
Nondimensional  time 

Nondimensional  solar  simulator  intensity,  — 
Nondimensional  volume 
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Equilibrium 

External 

Gas 

Convection  or  conduction  surface 

Internal 

Numerals 

Model 

Reference,  initial,  outside 
Prototype 

Reference,  radiation  surface 
Solar,  solid 
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SECTION  I 
INTRODUCTION 


Thermal  model  testing  of  space  vehicles  and  associated  components 
has  increased  because  full-scale  environmental  tests  have  become  more 
difficult  to  accommodate  in  existing  ground  test  facilities.  This  approach 
requires  that  a  thermal  model  must  be  designed  that  will  accurately  pre¬ 
dict  the  thermal  response  of  the  full-scale  prototype.  To  accomplish  the 
practical  design  of  thermal  models,  sufficient  experience  must  be  gained 
to  ensure  that  the  modeling  techniques  used  will  predict  the  thermal 
behavior  of  the  prototype. 

Thermal  model  testing  has  been  the  subject  of  several  papers  during 
the  last  few  years.  In  particular,  thermal  scale  modeling  was  discussed 
by  Vickers  (Ref.  1)  and  Wainwright  (Ref.  2).  A  thorough  treatment  of 
thermal  modeling  rules  was  given  by  Chao  and  Wedekind  (Ref.  3),  and 
Jones  (Ref.  4)  derived  an  extensive  set  of  similarity  ratios  to  be  pre¬ 
served  when  several  bodies,  each  having  its  own  characteristic  tempera¬ 
ture,  are  scaled  together. 

The  object  of  this  study  is  to  discuss  several  aspects  of  thermal 
model  testing  in  space  simulation  chambers.  The  subjects  discussed 
include:  a  criterion  for  the  application  of  isothermal  modeling  rules, 
the  application  of  nonisothermal  modeling  rules,  the  effect  of  internal 
convection  and  internal  radiation  on  the  thermal  balance  of  a  model,  and 
a  criterion  to  determine  the  degree  of  intensity  nonuniformity  for  a  solar 
simulator  in  terms  of  desired  temperature  accuracy.  By  treating  these 
aspects  of  thermal  modeling  separately,  and  correlating  this  with  previous 
work  (Refs.  1-4),  it  is  intended  that  a  practical  technique  of  thermal  model 
testing  will  be  implemented. 


SECTION  II 

SIMPLIFIED  THERMAL  MODELING 


The  heat  conduction  equation  in  body  number  i  of  an  assembly  of  n 
bodies  is 

Pi  c,  =  div  (Ai  Srad  Ti)  +  s.  (1) 

where  S;  is  the  internal  heat  generation  per  unit  volume.  Radiation  heat 
exchange  at  the  boundary  is  expressed  by  the  condition 


=  q 


i>  net 


b  o  u  n  d  a  ry 


(2) 
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where  Z-,  is  the  coordinate  normal  to  the  surface  and  q;  net  is  the  net 
radiative  heat  flux  per  unit  area.  The  general  form  of  this  heat  flux 
has  been  discussed  by  Chao  (Ref.  3),  The  quantity  q:  net  may  include 
contributions  from  solar  radiation,  albedo,  and  earth  radiation  as  well 
as  radiation  from  other  bodies  and  other  parts  from  the  same  body, 
externally  or  internally. 

In  this  report  the  external  prototype  and  model*  geometry  are 
assumed  to  be  similar.  Also  the  radiative  surface  properties  are  pre¬ 
served  in  scaling.  Furthermore,  if  the  temperature  at  homologous 
points  is  the  same  on  model  and  prototype,  effects  of  changes  in  surface 
emissivity  or  absorptivity  with  temperature  do  not  enter  into  the  scaling 
problem.  In  general,  this  state  can  be  achieved  only  if  the  model  has  the 
same  surface  finish  as  the  prototype. 

If  the  internal  heat  generation  term,  S;,  is  neglected,  the  non- 
dimensional  form  of  Eqs.  (1)  and  (2)  is 


1 

0 

-< 

! 

rS> 

7  <?3  d 3  L1 

I _ 1  _ _  1  - 

jL\J) 

dr  /Vc.L3  ^ 

Vt?^2  dr]3  dJ 

/ 

J 

(3) 

CTTa3d 

^  '  boundary 

Xj  qi 

] 

where 

T  =  0TO  X  =  £L 

t  =  n0  Y  =  77L 

q;  ^  qvT/  Z  =  <d 

The  coordinates  X  and  Y  are  along  the  surface  and  are  scaled  in  the  same 
ratio,  whereas  the  coordinate  Z  can  be  scaled  in  a  different  ratio,  allow¬ 
ing  a  minor  geometric  distortion. 

Vickers  (Ref.  1)  stated  that  it  is  not  possible  to  carry  out  thermal 
modeling  with  preservation  of  material,  surface  finish,  and  temperature 
if  the  model  has  complete  geometric  similarity.  This  is  shown  by  list¬ 
ing  the  similarity  parameters,  from  Eq.  (3),  to  be  preserved  between 
prototype  and  model: 


*The  term  prototype  is  used  to  designate  the  space  vehicle,  and  the 
term  model  refers  to  a  scaled  (larger  or  smaller)  version  of  the  space 
vehicle. 
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Since  it  is  extremely  attractive  to  preserve  both  temperature,  as 
explained  above,  and  materials  in  thermal  scaling,  minor  geometric 
distortions  are  introduced  where  possible.  If  the  temperature  field  is 
essentially  two-dimensional,  the  boundary  condition  in  Eq.  (3)  can  be 
introduced  in  the  differential  equation 


d01  A,t0 

(d2 

+  *'\6  + 

l3  ffT'-5 

3~  p  c .  L 2 

r  1  1 

ye 

»,*j  ■ 

,1  A  ^ 

1 

and  the  similarity  parameters  are 


or  for  material  and  temperature  preservation,  the  ratios 


have  to  be  preserved.  This  thin  shell  approximation  was  given  in  Refs.  2 
and  3  and  experimentally  verified  in  Ref.  5. 


If  the  bodies  or  parts  of  the  bodies  to  be  scaled  can  be  considered 
isothermal,  another  set  of  similarity  rules  can  be  derived,  affording 
more  freedom  than  given  by  Eq.  (5).  Isothermal  bodies  have  been  con¬ 
sidered  in  detail  by  Jones  (Ref.  4).  Here  a  simple  derivation  is  given, 
starting  from  Eq.  (1)  and  assuming  at  the  onset  that  external  radiation 
scaling  is  satisfied  by  using  similar  geometries  and  identical  surface 
properties  and  neglecting  St.  The  right-hand  side  of  Eq.  (1)  now  con¬ 
sists  of  two  parts  -  (1)  heat  transfer  through  a  surface  Ah  from  another 
body  j  adjacent  to  the  body  i  with  a  heat  transfer  coefficient  H;j  and 
(2)  radiation  transfer  through  a  surface  Ar  with  the  flux  q,  =  q;_  net 
replacing  the  boundary  condition  Eq.  (2).  Therefore,  for  n  bodies. 


d  T; 


j=n 


A  ",  (T,  -  T.)  -  „  A„ 


1*  1 


Using  the  nondimensional  notation  as  before  and  in  addition 


V-  =  <5,  Vj 

Ahi 

=  abi  ahi 

Hu  =  Vu 

“  °r.  ar. 

Equation  (6)  becomes 

d0,  /  hJ;  ahl 

r  h  cii 

ft  -  *«) 

i  dr  1  p.  v 

,  )  i=i  "  hl 

[  P,c, 

)Ai 

%  ar 


(6) 


(7) 
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and  the  similarity  parameters  to  be  preserved  are 


ahi  lo 

p,  C  V 
r  l  11 


and 


1°  Tt/  ar i 

\  P,C,V, 


Material  and  temperature  preservation  reduce  this  further  to 

(^) 

If  the  prototype  is  scaled  by  a  factor  K,  the  external  surface  area,  ar;, 
will  be  scaled  by  a  factor  K2  for  geometric  similarity  required  by 
radiation  exchange.  The  last  parameter  in  Eq.  (8)  gives,  then, 

_  (°“)p  __  1  (V^ro 

(*-)P  (ari)m  <vi)p  (vi)p 

This  introduces  some  freedom  in  selecting  the  time  scale: 


(a)  If  (t„)m  =  (tQ)p  is  selected,  the  volume  must  be  scaled  by  K2, 

which  is  possible  by  leaving  the  material  thickness  unchanged. 
This  was  done  in  Ref.  6  where  the  volume  (and  heat  capacity) 
was  scaled  by  a  factor  K2.  The  first  parameter  in  Eq.  (8) 
gives,  now, 

(hjj  ahi)m  -  K  (hjj  ahi)p 

In  general,  one  dimension  of  the  cross  section,  ahi»  will  be 
along  the  external  surface  and  one  dimension  will  be  in  the 
direction  of  the  thickness,  so  that  (ah , ) ^  =  KUhi)  and  therefore 

(hn>„  -  *  0>i.)p 

The  model  has  to  be  divided  into  isothermal  parts  separated  by 
junctions  with  heat  transfer  coefficients  smaller  than  the 
original  ones  in  the  case  of  K  <  l;  i.  e.,  additional  thermal 
resistance  has  to  be  introduced.  When  there  is  only  radiative 
heat  transfer  between  the  isothermal  parts,  only  the  last 
parameter  in  Eq.  (8)  has  to  be  satisfied. 


(b)  If  (t0)m  =  K(t0)p  is  selected,  the  volume  scaling  is  (Vj)m  =  K3(Vj)( 
and  the  model  is  scaled  equally  in  all  directions.  The  heat 
transfer  coefficients  are  now  scaled  according  to 


(hij)  (v,)  (l0)  (ahi)  K  K 

p  P  m  jh 

The  heat  transfer  coefficients  remain  unchanged. 
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(c) 


If  (to)m  =  KJ(t0)p,  the  volume  is  scaled  by  K*  and  the 
heat  transfer  coefficients  are  scaled  by 


(hi;) 


m 


(fcij) 


P 


(vi)  <‘o>  (ah.)  „  .  , 

_ m  _ P_  _ -  K  ■  —  .  — 

(vi)  (h.)  (ah0  K’  K> 

pm  m 


1 

K 


This  is  the  same  result  as  that  obtained  with  the  thin  shell 
approximation  where  the  thermal  resistance  is  scaled  by 
a  factor  K  automatically. 


SECTION  III 
ISOTHERMAL  BODIES 


In  order  to  apply  isothermal  modeling  rules,  it  has  to  be  decided 
when  a  body  or  part  of  a  body  can  be  considered  isothermal.  If  at  any 

time  the  ratio  4j?-  «  1 ,  the  body 


is  said  to  be  isothermal.  Here 
AT  is  a  typical  temperature  differ¬ 
ence,  say,  between  the  hottest  and 
coldest  point  on  the  body.  An 
equivalent  statement  is 


dTi  i 

^i«l 


(9) 


where  Li  is  a  typical  length  in  the  direction  of  X, .  It  follows  from  the 
heat  conduction  Eq.  (1)  that 


T  -  0  (x-  L  If)' 


(10) 


It  is  assumed  that  the  internal  heat  generation  can  be  neglected.  Com¬ 
bining  Eqs.  (9)  and  (10)  leads  to  the  condition  to  be  satisfied  for  the 
isothermal  assumption: 


l  dr 
T  di 


«  1 


or,  also. 


«  1 


pc  L 
A  h 

where  t0  is  a  typical  time,  for  instance  the  period  in  a  test  with  cyclic 
solar  simulation. 


(11) 


*The  notation  0  ( . )  means:  on  the  order  of  ( 
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Condition  (11)  applies  to  the  overall  length  of  a  body  shown  in  the 
accompanying  sketch.  This  condition  is  certainly  satisfied  for  a  minor 

dimension,  such  as  the  thickness 


<3.  However,  if  solar  radiation 
in  one  side  of  the  body  provides 
.net  the  heat  input,  the  boundary  condi¬ 
tion  of  Eq.  (2)  gives  an  additional 
requirement.  From  Eq.  (9), 
gfad  (T)  •  d  <<;  j 


f  q.  )  d. 
*  net'  . 

XT 


«  1 


(12) 


As  an  example,  the  requirements  for  isothermal  bodies  are  listed 
below,  with  T  =  300°K,  nel  =  0.14  watt/cm2  (solar  radiation),  and  t0  =  100  min 
(orbital  period): 


Material 

Copper 

Aluminum 

Stainless  Steel 

X,  watt/cm  -°K 

4.  0 

1.  0 

0.  13 

A/p  c,  cm2/ sec 

1.  1 

0.44 

0.  04 

d,  cm  « 

8600 

2150 

280 

L  ,  cm  « 

6600 

2640 

240 

It  is  concluded  from  this  table  that,  even  for  materials  with  a  high 
thermal  conductivity,  truly  isothermal  conditions  are  only  achieved  for 
very  small  bodies,  say,  with  a  length  L  =  10  to  30  cm  for  copper. 

The  isothermal  criterion  is  illustrated  by  experiments  reported  in 
Ref.  6.  Three  bodies  -  a  plate,  a  cylinder,  and  a  sphere  -  were  placed 
in  a  cold  vacuum  chamber  (Fig.  1).  In  this  case  the  plate  was  heated 
periodically.  The  extreme  temperatures  of  the  aluminum  cylinder, 
heated  by  radiation  from  the  plate,  were  equal  within  a  few  degrees 
(Fig.  la).  On  the  other  hand,  the  extreme  temperatures  of  the  stainless 

steel  sphere  differed  as  much  as  140°K  (Fig.  lb).  The  parameter  — 

A 

was  equal  to  0.  06  in  the  former  and  1.  6  in  the  latter  case.  Since  in  the 
latter  case  the  isothermal  condition  was  not  met,  a  scaling  procedure 
based  on  the  assumption  that  each  of  the  three  bodies  was  isothermal 
did  not  result  in  an  accurate  prediction  of  prototype  temperatures  from 
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a  thermal  model  experiment.  This  point  was  demonstrated  in  Ref.  6. 
Even  an  aluminum  sphere  would  not  quite  meet  the  isothermal  require¬ 
ment  [Eq.  (11)]  .  * 


Scale,  cm 

I _ I _ l _ I _ l 

0  20  40 


Aluminum  Cylinder  — ^  Plate  Heated  Periodically  with 

Internal  Heater,  Ref.  6 


Prototype  Data  Only  Presented 


a.  Aluminum  Cylinder 


b.  Stainless  Steel  Sphere 

Fig.  1  Experimental  Verification  of  Isothermal  Criterion 


*In  Ref.  6,  test  results  with  an  aluminum  sphere  are  given,  but  in 
that  test  the  sphere  also  had  an  internal  heater  to  help  produce  uniform 
temperatures. 
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SECTION  IV 

NONISOTHERMAL  BODIES 


When  the  same  configuration,  with  a  stainless  steel  sphere,  was 
scaled  using  the  thin  shell  approximation,  i.  e. ,  scaling  was  carried  out 
on  the  basis  of  thermal  resistance  rather  than  thermal  capacity,  good 
agreement  between  homologous  temperatures  on  the  model  and  proto¬ 
type  was  obtained  (Ref,  5). 

Another  example  of  the  use  of  the  thin  shell  (two  dimensional)  model¬ 
ing  rules  given  by  Eq.  (5)  is  shown  in  Fig.  2,  where  typical  temperatures 
of  a  prototype  and  a  half-scale  model  are  compared.  In  this  case,  the 
bodies  were  heated  periodically  with  a  carbon  arc  solar  simulator  in  a 
cold  (77°K)  vacuum  chamber. 
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SECTION  V 

INTERNAL  CONVECTION 


In  thermal  testing  of  space  vehicles,  the  gas  (or  liquid)  in  the 
vehicle  may  significantly  influence  the  heat  balance  and  the  temperature  - 
time  history.  Tests  carried  out  in  a  space  chamber  with  gravity- 
induced  convection  do  not  simulate  free-space  conditions.  Only  a  facility 
for  testing  free-fall  models  would  properly  simulate  zero  gravity  condi¬ 
tions.  While  it  is  clear  that  the  liquid  in  a  thin-walled  tank  influences, 
markedly,  the  temperature-time  history  of  such  a  tank  when  it  is  sub¬ 
jected  to  a  periodic  solar  source  in  cold  space,  it  is  not  immediately 
apparent  when  the  heat  capacity  of  a  gas  and  the  heat  transmission 
through  the  gas  are  important.  The  following  derivation  is  therefore 
given  primarily  to  establish  the  limits  for  negligible  convection  effects. 

The  heat  balance  equations  for  the  gas  and  a  solid  container  of 
arbitrary  shape,  as  sketched  below,  are 

<M„A,  -  t  o  Ta*  Ae  -  h  (T,  -  Tg)A,  -  (13) 


(T-  -  Te) 


A.  = 


TT 


(14) 


The  notation  is  shown  in 
the  sketch.  The  solid  is 
considered  isothermal 
(infinite  thermal  conduc¬ 
tivity).  The  coefficient  h 
applies  to  the  overall  heat 
transfer  between  the  solid 
and  the  gas.  The  container 
radiates  to  cold  (0oK)  space. 
The  equilibrium  tempera¬ 
ture  is  given  by 

a  A,  =  e  cr  Tr-4  A„ 


where  As  is  the  projected  surface  area  normal  to  the  solar  beam 
direction  and  Te  is  the  final  equilibrium  temperature. 
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Equations  (13)  and  (14)  can  be  written  in  nondimensional  form: 


P  ((1  +  Q)4  -  (dB  +  Q)4}  -  (0S  -  0g)  = 


(15) 


where 


and 


C 


c„  M„ 


es  ~  e g  =  c 


p  = 


f  cr  Ac 
h  A; 


<Te  -  T0)5 


(16) 


(17) 


This  set  of  equations  was  solved  numerically  for  the  following  boundary 
conditions: 


At  r  =  0  :  6B  =  Bs  =  0 
At  r-»oc  :  =  0S  =  } 

A  typical  solution  is  shown  in  Fig.  3.  The  results  of  calculations  for 
various  values  of  C,  P,  and  Q  are  summarized  in  Fig.  4,  which  is  a 
plot  of  the  decrease  in  at  0S  -  0.8  that  was  attributable  to  the  pres¬ 
ence  of  the  gas  versus  C  (the  ratio  of  thermal  capacity  of  the  gas  to 
that  of  the  container)  with  P  and  Q  as  parameters. 

For  a  sphere  of  1-m  diameter  filled  with  air  at  1  to  100  atm  at 
room  temperature,  it  was  estimated  that  h  =  3  to  30  watt/m3  °K.  With 
<  =  1.0  and  Ae/Ai  =  1,  this  gives  P  =  0.05  to  0.5.  A  metal  sphere  of  1-m 
diameter  and  thickness  of  one  hundredth  of  its  diameter,  filled  with 
air  at  1  atm  has  a  value  of  C  on  the  order  of  0.  01. 

For  the  cases  considered,  A 0S  <  0.5C,  and  thus  for  "ordinary"  con¬ 
ditions  it  is  expected  that  the  effect  of  the  gas  on  the  temperature  of  the 
solid  is  very  small  (a  few  degrees)  and  probably  within  the  measuring 
accuracy  for  many  tests.  Under  space  conditions,  the  heat  transfer 
coefficient,  h,  is  expected  to  be  considerably  smaller  than  in  a  space 
chamber;  therefore,  P  will  be  larger,  which  tends  to  reduce  the  effect 
of  the  gas  further. 

Figure  4  helps  to  set  the  limits  when  convection  can  be  neglected. 
However,  it  is  interesting  to  consider  convection  scaling  parameters. 
From  Eqs.  (15)  and  (16)  it  follows  that  C,  P,  and  Q,  together  with  the 
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Fig.  4  Effect  of  Gas  on  the  Temperature  of  a  Solid  during  Transient  Conditi 


AEDC-T  R-67-38 


c  M 

characteristic  time  tr  =  g  '  8  are  the  similarity  parameters.  Making 

h  A; 

the  same  restrictions  as  before  -  temperature  and  surface  properties 
preservation  -  the  similarity  parameters  become 

A«  \  AA 

h  A;  /  ’  [t.h  A, 


Since  for  similarity  hAj  -  Ae,  one  can  also  write 


(18) 


The  first  parameter  appears  to  cause  no  special  problems.  Since  it  is 
expected  that  under  space  conditions  h  is  much  smaller  than  under  labora¬ 
tory  conditions,  the  internal  area  of  the  model  has  to  be  reduced  or  h  has 
to  be  reduced  by  placing  internal  low  conductivity  partitions  in  the  model. 
If  in  the  third  parameter  cs  remains  unchanged,  then  for  a  thin  shell 
Ms  -  K1  (surface)  x  KJ  (thin  shell)  and  A4  -  K2  so  that  the  time  scale  is  propor¬ 
tional  to  K2.  On  the  other  hand,  if  the  body  were  scaled  isothermally, 

Ms  -  K2  and  the  time  would  not  be  scaled.  Thus,  convection  scaling  does 
not  lead  to  requirements  conflicting  with  other  thermal  scaling  param¬ 
eters,  at  least  under  the  present  simplifying  assumptions. 


SECTION  VI 
INTERNAL  RADIATION 


To  assess  the  importance  of  internal  radiation  in  thermal  testing, 
the  temperature  distribution  of  a  sphere  with  an  internally  radiating 
surface  is  studied.  When  a  hollow  sphere  is  subjected  to  solar  radia¬ 
tion,  as  sketched  in  Fig.  5,  the  temperature  distribution  is  given  by 
the  following  differential  equation: 


<?2T  dT  (iB1  (7RJ  4  f.ffR3  f  4 

-  +  cot  v -  +  -  l(v)  qD  cos  V  -  +  (.1 - T  +  — -  I  T  sm 

du 3  dv  Ad  8  1  *  Ad  2Ad  Jo 


R  dT 
fdf  -  pc  —  — - 
A  a  t 


0 


By  introducing  the  nondimens ional  variables. 


6  =  —  where  T(l) 

T  (1) 


temperature  at  v  =  0 


T 


~~  »  with  t0 

^  n 


p  cRJ 
A 
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Hollow  Sphere 


a  -  Absorptivity 

cos  v  a  M 


6  = 


AdT(l) 


T(l) 


1/4 


4  +  6  /  ct  \  /  qs\ 
4  +  46  Vig/  \a~J 


4— 


Fig,  5  Steady-State  Temperature  Distribution  of  a  Sphere  with  Internal  Radiation 
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and  using  the  notation 


p  =  COS  V 


this  reduces  to 


A.  <1  T(l) 

(  l  +  S/}  a  q  R} 

‘9 


y 


(1 


n  +  h 
]  +  <v4 


(19) 


The  function  f(fd  is  f(jx)  =  1  for  0  <  p  <  1  and  lip)  =  0for-l<pi<0  ,  and 
accounts  for  the  directional  effect  of  the  solar  source.  The  derivation 
of  Eq.  (19)  is  rather  long  and  is  presented  here  only  to  bring  out  the 
influence  of  internal  radiation. 


For  thermal  scaling,  the  parameters  y  and  8  have  to  be  preserved. 

If  surface  properties  are  preserved,  the  reference  temperature  T(l)  is 
preserved;  and  if  material  properties  are  preserved,  the  parameters  y 

.  n 

and  8  reduce  to  the  requirement  that  be  preserved.  Thus,  if  =  K, 

R  Rp 

the  thickness  is  scaled  by  *  K2  and  the  time  is  scaled  by 

/  v»  dP 

(^m  =  /Hrl\  _ 

(lo)p  \RP/ 

This  is  the  same  scaling  result  as  that  obtained  with  the  thin  shell 
approximation,  since  the  same  assumptions  are  involved.  Without 
going  into  detailed  numerical  solutions,  it  suffices  here  to  plot  the 
temperature  distribution  along  the  surface  of  the  sphere  for  the  two 

cases;  8  =  =  0  and  1  (Fig.  5).  The  case  of  y  =  0  corresponds  to  zero 

^  e 

conductivity  in  the  sphere  wall,  and  y  =  «■  is  the  case  of  infinite  conduc¬ 
tivity  or  an  isothermal  sphere.  It  is  concluded  that  the  internal  radia¬ 
tion  properties  of  a  spacecraft  model  would  have  to  be  reproduced 
accurately  in  model  testing,  in  particular  where  complicated  internal 
arrangements  occur  since  corrections  would  then  be  difficult  to  make. 


SECTION  VII 

SOLAR  SIMULATOR  UNIFORMITY  CRITERION 


The  degree  of  simulation  achieved  by  a  solar  simulator  is  a  func¬ 
tion  of  how  accurately  certain  characteristics  of  the  actual  solar  source 
are  simulated.  One  of  these  characteristics  is  the  uniformity  of  the 
solar  simulator  beam  intensity. 
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In  a  general  way,  the  degree  of  nonuniformity  that  can  be  tolerated 
in  a  simulated  solar  beam  depends  on  the  test  model  geometry,  the 
model  material,  and  the  allowable  error  in  the  temperature  distribu¬ 
tion.  A  simplified  case  v/ill  be  considered  as  an  initial  step  in  the 
development  of  a  general  solar  simulator  intensity  uniformity  criterion. 


Since  nonuniformity  effects  are  expected  to  be  most  severe  for  thin 
shells  and  since  most  spacecraft  are  constructed  of  plate  material, 
consider  the  heat  balance  for  a  flat  plate  exposed  to  radiation  in  the  z- 
direction.  Assume  no  heat  flow  in  the  y-  and  z-  directions.  The  former 
implies  a  two-dimensional  geometry  of  both  the  model  and  the  solar 
intensity  nonuniformities,  whereas  the  latter  is  the  thin  shell  approxi¬ 
mation.  With  these  assumptions, 

pc  11  =  A  +  -1  (q  -  tof) 

<?t  ax2  d  1 


z 


where  d  is  the  material  thickness.  Using  the 
nondime  ns  ional  variables  introduced  in  Eq.  (3) 
for  T,  t,  and  X,  while  introducing  the  addi¬ 
tional  nondimens  ional  variable  q  =  <f>  qs> 
gives 


66  _  A  tD\  6*6 
dr  \pc  If  d^1 


ea  T  st 

_ O  0 

p  c  d 


6 4 


(20) 


As  a  first  approximation  to  a  general  nonuniformity  criterion  for 
solar  simulation,  Eq.  (20)  will  be  solved  in  a  simplified  form.  For 
this,  assume  that 


(1)  AA.  -  o,  steady-state  condition, 

dr 

(2)  the  last  term  in  Eq,  (2)  is  negligible,  and 

(3)  d>  =  <M£) .  In  general,  6  =  6  ( £,  r) . 


With  the  above  assumptions,  Eq.  (20)  becomes 

d'd  q  iS  ,  4 

—  +  r  <p  =  o  ;  r  =  —  (2i) 

df  AdT0 

Since  <£  is  equal  to  the  ratio  of  the  simulated  radiation  intensity  to 
the  average  solar  radiation  constant,  the  term  F  can  be  called  the  intensity 
nonuniformity  parameter.  The  magnitude  of  this  term  dictates  the  non¬ 
uniformities  in  the  solar  simulator  intensity  that  can  be  tolerated  for  given 
allowable  temperature  tolerances. 
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As  an  illustration,  assume  that  4>  =  (1  +  A  sin  £??).  The  degree  of 

nonuniformity,  A  =  -  1 ,  is  as  shown  in  Fig.  6.  Equation  (21) 

becomes 

-d  ^  +  r  (1  +  A  sin  £  rr)  =  0 


1 

ISl 
c r 

O" 

n 

0 


0 

1 


0 


Fig.  6  Solar  Simulator  Nonuniformity 

Considering  the  case  of  a  plate  with  a  constant  temperature  at  both  ends 
(Fig.  6),  the  solution  is 

d  =  AL  sin  ^  _  _L  r  £ 1  +  -2-  r  |  +  i 

n2  2  2 

The  dashed  line  on  the  graph  in  Fig.  6  is  the  ideal  temperature  distribu¬ 
tion  with  a  uniform  beam  intensity  (A=0)  and  the  solid  line  is  the  distribu¬ 
tion  with  a  nonuniform  intensity  (A*0,  n^O).  The  error  in  temperature 
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at  any  particular  point  is 

E  =  (0Auo  -  (0)A  =  O  =  sin  £* 

and 


(22) 


Since  n± = 


L,  it  follows  from  Eqs.  (21)  and  (22)  that 


E 


max. 


AT  1  q  L2 
T  o  »AdT0 


(23) 


As  n  approaches  infinity,  Eraax.  approaches  zero  regardless  of  A.  This 
shows  that  if  the  model  dimensions  are  large  compared  to  the  size  of 
the  intensity  nonuniformities,  the  error  in  temperature  becomes  small. 
Assuming 


L  =  100  cm,  n  =  10,  S.  =  10  cm,  d  =  0.3  cm,  and  qs  =  0.14  watt/ cm2 

gives 

q  /r)2n2d\ 

A  =  ma*-  -  1  =  - H  A  A  T  =  0.221  AAT 

%  \%h  / 

So  that  per  degree  Kelvin  error  (AT  =  PK),  the  maximum  allowable 
intensity  nonuniformity  is  as  shown  below; 


Copper 

Aluminum 

Stainless  Steel 

A,  watt/ cm— °K 

4.0 

1.  0 

0.  13 

A/AT,  percent 

±88.  4 

±22.  1 

±2.  9 

The  importance  of  geometry  of  the  model,  the  solar  intensity  uni¬ 
formity,  and  the  material  properties  are  demonstrated  by  Eq.  (23)  and 
the  table  above. 

Thermal  modeling  tests  have  been  conducted  with  a  carbon  arc 
solar  simulator  with  two  different  model  configurations.  The  same 
solar  simulator  was  used  in  both  cases.  The  models  in  one  test  were 
fabricated  from  stainless  steel,  whereas  the  models  used  in  the  second 
test  were  fabricated  from  aluminum  sheet.  The  test  articles  fabricated 
with  aluminum  had  a  thickness  that  was  double  the  thickness  of  the 
stainless  steel  models,  resulting  in  a  higher  value  of  thermal  conduct¬ 
ance.  Referring  to  the  intensity  nonuniformity  parameter  as  incorpo¬ 
rated  in  Eq.  (23)  and  considering  that  the  same  solar  simulator  was 
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used  on  both  series  of  tests,  gives 

(AT)s.s.  (Xd)s,s,  =  0.0)  (0.125)  =  Ug 

(AT)a1  (Ad)Al  (0.13)  (0.06S) 

This  shows  that  the  error  that  could  be  expected  between  analogous 
temperatures  on  the  scale  model  and  prototype  would  be  approximately 
15  times  greater  with  the  test  articles  made  of  stainless  steel  than 
for  those  fabricated  from  aluminum.  Typical  temperatures  are 
plotted  in  Fig.  7  for  the  stainless  steel  and  aluminum  test  articles, 
where  this  ratio  is  only  five.  It  is  noted  that  the  last  term  in  Eq.  (20) 
tends  to  reduce  the  effects  of  intensity  nonuniformities  on  temperature 
distribution.  The  effect  of  this  term  and  two-dimensional  heat  con¬ 
duction  during  time-periodic  heating  contribute  to  the  above  discrepancy. 


SECTION  VIII 
CONCLUDING  REMARKS 


The  emphasis  in  this  report  has  been  on  thermal  modeling  with 
simultaneous  preservation  of  temperature,  surface  properties,  and 
materials.  It  is  felt  that  for  a  complicated  model  this  is  the  most 
practical  method.  When  the  bodies  cannot  be  divided  into  isothermal 
parts  in  a  convenient  manner,  the  nonisothermal  scaling  may  introduce 
impractical  material  thicknesses,  particularly  when  the  model  is  to  be 
significantly  smaller  than  the  prototype.  In  such  a  case,  much  can  be 
learned  from  partial  scaling,  for  instance  on  two  or  three  different 
models,  and  extrapolation  to  the  correct  scaled  size.  Tests  with  two 
or  three  different  materials  may  also  provide  useful  data  for  this 
extrapolation. 
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Prototype 


0  120  240 

tp,  min 


Fig.  7  Experimental  Data  Illustrating  Solar  Simulator  Nonuniformity  Effects 
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